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Against expectations, robust switchable ferroelectricity has been recently observed in ultrathin 
(1 nm) ferroelectric films exposed to air [V. Garcia et al., Nature (London) 460, 81 (2009)]. Based 
on first-principles calculations, we show that the system does not polarize unless charged defects or 
adsorbates form at the surface. We propose electrochemical processes as the most likely origin of 
this charge. The ferroelectric polarization of the film adapts to the external ionic charge generated 
on its surface by redox processes when poling the film. This, in turn, alters the band alignment at 
the bottom electrode interface, explaining the observed tunneling electroresistance. Our conclusions 
are supported by energetics calculated for varied electrochemical scenarios. 



I. INTRODUCTION 

Complex oxides have long been viewed as possible 
candidates for the next generation of electronic devices, 
which require reduced feature sizes, enhanced operating 
speeds and low consumption. Amongst oxides, ferro- 
electrics offer the ability to store information in a non- 
volatile manner via fast reversible polarization switching 
in ferroelectric random-access memory (FeRAM). The 
observation of giant tunneling electroresistance (TER)^ 
in ultrathin (3 unit cells) ferroelectric films has recently 
opened a novel paradigm for device design based on these 
materials^ ^. 

Although the experiments^ ascribed TER to ferroelec- 
tricity, which appeared robust and switchable, how the 
polar state is stabilized in such thin films is by no means 
established. In principle, a ferroelectric film with an ex- 
posed surface cannot sustain a monodomain polarization 
perpendicular to the surface, because of the strong de- 
polarizing field that would inevitably arise^^. Charged 
particles from the environment could in principle cancel 
the depolarizing field"^ (Fig. 1 left). So far, however, the 
only chemical control of switching in air relates to neutral 
species, 02^^^ (Fig. 1 center). It is then not clear how 
neutral gas-phase molecules could interact with a biased 
atomic force microscopy (AFM) tip to produce the polar 
state. 

Here we argue that the voltage applied with the AFM 
tip induces electrochemical switching (Fig. 1 right), i.e. 
redox processes that are essential to liberate free charge 
and therefore screen the depolarizing field. This process 
would act as a nanobattery, rather than a nanocapaci- 
tor. Note that the same mechanism could explain other 
effects at oxide interfaces, such as the switchable two- 
dimensional electron gas (2DEG) at the LaAlOa/SrTiOa 
interfac^^, where the switching appears to be mediated 



by surface chargd^. 

To explore this mechanism we consider the system 
studied experimentally in Ref. 1, consisting of a com- 
pressively strained nanometer-thick BaTiOs (BTO) film 
on a Lao.rSro.sMnOs (LSMO) bottom electrode. Here 
we show, using first principles calculations, that (i) the 
pristine system (clean BTO surface with an ideal Ti02 
termination) does not allow for a ferroelectric polariza- 
tion, P, normal to the surface despite the large compres- 
sive strain; (ii) a non-zero P is crucially dependent on 
the presence of a surface external ionic charge, in the 
form of defects or adsorbates; and (iii) the energetics 
for the formation of oxidized or reduced surface defects 
support the electrochemical switching model. We also 
find (iv) a systematic change in band offset with screen- 
ing charge density, which we identify as the microscopic 
mechanism behind the experimentally observed TERP, 
and (v) a large magnetoelectric coupling, due to the ac- 
cumulation or depletion of spin-polarized carriers at the 
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FIG. 1: Schematic illustration of the conventional (left) and 
redox (center) mechanisms for ferroelectric screening in the 
absence of a top electrode. The presence of a biased tip can 
promote an alternative redox mechanism that provides an ex- 
ternal circuit for the screening electrons (right). 
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interface with ferromagnetic LSMO. The connection be- 
tween these effects can be summarized as fohows: under 
open-circuit boundary conditions the electric displace- 
ment field D within the film, the change in magnetization 
at the interface AM and the interface dipole, are all pro- 
portional (or equal) to the external ionic charge density, 
Q per unit surface produced by the redox processes. 



II. METHODS 

The density-functional theory (DFT) calculations are 
performed using the spin-polarized Wu-Cohen (WC) 
exchange-correlation functional^ ^, as implemented in the 
Siesta cod ^^^ [ We find GGA-WC to reproduce 
buUpl and surfacd^^ properties of LSMO that were 
calculated using the Perdew-Burke-Ernzerhof (PBE) 
scheme^^; at the same time, GGA-WC is more appro- 
priate for ferroelectric oxides. The LSMO/BTO system 
consists of 5.5 unit cells of LSMO (Mn02-terminated) 
stacked with 3 unit cells of BTO along the c direction 
in a slab geometry. The supercell contains a 15 A thick 
vacuum layer and has either 2x2 or V^xV^ in-plane 
periodicity (see Fig. 2). The 5.5 unit cells of LSMO are 
thick enough to show bulk-like features in the center, and 
3 unit cells of BTO was experimentally shown to be thick 
enough for ferroelectricit}'^ . We use a dipole correction 
to simulate open-circuit boundary conditions, enforcing 
zero macroscopic electric field in the vacuum layer. We 
constrain the in-plane lattice parameter to experimen- 
tal bulk NdGaOs (NGO) to reproduce the experimental 
conditions of Ref.^; this imposes a large (3%) compressive 
strain on BTO. Based on this slab geometry, we perform 
a number of calculations where we vary the surface com- 
position by introducing defects or adsorbates. In partic- 
ular, we simulate the clean Ti02-terminated surface (we 
shall refer to this structure as "pristine" henceforth) ; one 
O vacancy ("0-vac") or adatom ("0-ads") per 2x2 sur- 
face ceh; one Hadatom ("H") or OH group ("OH") per 
a/2 X ceh Hereafter we shall discuss the results 
with special regard for the presence or absence of ferro- 
electric polarization in each case. 



III. DISCUSSION 

A. The pristine system 

Fig. 2 shows the relaxed out-of-plane structural dis- 
tortions as a function of the surface chemical environ- 
ment. The pristine system is characterized by negligi- 
ble distortions in the interior of the BTO film, suggest- 
ing the absence of macroscopic P in this system. Only 
a surface rumpling is present, resulting in a small net 
inwards dipole (non-switchable) that decays rapidly to- 
wards the bulk (a surface rumpling is a known general 
feature of oxide surfaces, in particular the Ti02 termina- 
tion of BTO^^ ^^). A vanishing P is consistent with the 
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FIG. 2: (Color online) Cation-anion splittings, 5z — 2;cation — 
^anion through the LSMO/BTO slab (the bottom half of 
LSMO is not shown). The dotted lines correspond to the 
average of the AO and BO2 layer anion-cation splitting for 
the inwards and outwards P in bulk BTO, strained to NGO. 



open-circuit boundary conditions, despite the large com- 
pressive strain. In absence of a top electrode the macro- 
scopic electric displacement field D in BTO is equal and 
opposite to the density of external surface charge. As 
this charge is zero at the clean Ti02 surface, the film is 
constrained to a paraelectric state. 



B. Chemical switching 

To illustrate possible screening scenarios, we now in- 
clude representative surface defects^^. The O-vac and 
O-ads systems are both characterized by large ferroelec- 
tric distortions (Fig. 2). These are comparable to the 
strained bulk, where we calculate a spontaneous polar- 
ization Po=0.369 C/m^ (0.35 e/S). This result is again 
consistent with the constraint that D = —Q / S . In fact, 
one oxygen defect for every 2x2 unit cells (0.5 e/S) yields 
a larger surface charge than what would be sufficient to 
screen Pq- This justifies the larger cation-anion rum- 
plings that we obtain in the film compared with the bulk 
(Fig. 2). OH and H adatoms (with >/2x>/2 coverage to 
maintain Q / S) produce distortions of similar magnitude 
(Fig. 2). This confirms the generality of the ferroelectric 
switching mechanism: the ferroelectric state really de- 
pends on the net surface charge, and not on the chemical 
identity of the adsorbed species. 

In order to study the electrochemical switching (Fig. 
1 right), we commence by analyzing chemical switching 
(Fig. 1 center). Both are controlled by redox processes 
that transform bound charge into free charge, allowing 
for an electronic transfer between the surface defect and 
the metal substrate, but have different associated chem- 
ical sources/drains and energetics. Chemical switching 
was recently shown in a system consisting of PbTiOa on 
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FIG. 3: (Color online) Change in Mn magnetic moment (left), 
AM, near the interface (the interface Mn02 layer is on the 
right). Open symbols represent the change in 3d Cg occu- 
pation, and closed symbols the total magnetic moment. A 
schematic illustration (right) of the effect of 0-vac and 0-ads 
on the BTO polarization (arrows) and the Mn 3d Cg occupa- 
tion (blue lobes) and total Mn magnetic moment (numbers). 



SrRuO20^ and BTO films on Au or vacuunf. 

To assess whether these redox reactions are ther- 
mo dynamically accessible in typical experimental con- 
ditions, we estimate the formation energy of the de- 
fective systems taking the reactions: 1) Slab (pristine) 
Slab(0-vac)+l/202 and equivalent for O-ads, 2) 
l/2H20+l/402+Slab(pristine) ^ Slab(OH) and 3) 
1/2H20+Slab(pristine) ^ Slab(H)+l/402. The chemi- 
cal potential of the relevant molecular species is set to the 
calculated total energy of the spin-polarized molecule in 
a large cubic box. The results are summarized in Table 
1. They suggests that, whilst the oxygen adatom is likely 
to form under oxygen-rich conditions, the formation en- 
ergy for the oxygen vacancy is possibly too high to form 
even in oxygen-poor conditions. The calculated OH and 
H formation energies suggest that water is a very likely 
redox intermediate. Note that H2O is ubiquitous in most 
experiments performed in air, and was recently found to 
play a crucial role in AFM experiments performed on 
LaAlOa/SrTiOa^. Since both sets of reactions involve 
oxygen, we therefore expect that altering the surround- 
ing oxygen partial pressure would affect the stability of 
reduction or oxidation processes, consistent with the re- 
cently observed chemical switching^. 



TABLE I: The formation energy, Ef^ of the defective systems 
for O2 and H2O rich conditions (see text for definitions). 
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C. Electrochemical switching 

Now we discuss how the electrochemical processes 
could proceed in practice during the AFM switching ex- 
periments of Ref. |T] (general electrochemical processes on 
oxide surfaces are reviewed in Ref. 20 21). As schemat- 
ically shown in Fig. 1 (right), a biased tip close to con- 
tact can remove surface ions. These would then undergo 
a redox reaction at the tip surface. This process is fa- 
vored by the energy associated with the biased external 
circuit, QVexti but costs an energy equal to the change 
in binding energy of the ion to the ferroelectric surface 
and to the tip surface, A£^j (this effectively redefines the 
relevant chemical potential). By minimizing the Gibbs 
free energy of the system (see e.g. Ref. [22| or |23|) it 
can be shown that poling can stabilize redox defects if 
Vext > ^Ef/Q^ after which the equilibrium redox charge 
density, Q/S^ and polarization both grow with Vext- This 
electrochemical process would then act as a nanobattery, 
rather than a nanocapacitor. By controlling the envi- 
ronment (species and chemical potential) and Vext^ one 
may be also able to selectively control the active redox 
reaction, potentially opening new routes to surface redox 
catalysis. After removal of the tip, the surface redox den- 
sity from poling can remain, since the reverse reaction is 
now blocked by key reactants being removed with the 
tip. This would explain the observation of Ref. W that 
the domains are stable for a very long time after "writ- 
ing". Of course, lateral charge diffusion across domain 
boundaries^ may still occur in principle, but kinetic bar- 
riers are likely to hinder such processes. Therefore the 
bulk polarization, Pq, is expected to be an estimate of 
the equilibrium polarization after poling. We note that 
unlike in the LaAlOs/SrTiOs system where the polar- 
ization is driving the surface chemistr};^, in ferroelectric 
films we expect it is the surface chemistry (and poling) 
that is driving the polarization. This is because the en- 
ergy scale for changing the polarization is much larger in 
LaAlOs than in the ferroelectric. 



D. Magnetoelectric coupling 

The electronic transfer mechanism can be quantitively 
estimated through the change in magnetization of LSMO. 
LSMO is a half-metal with only Mn 3d Cg majority spin 
levels around the Fermi level. As the screening carri- 
ers are fully spin-polarized, an electronic transfer be- 
tween LMSO and the BTO surface results in a system- 
atic change of the magnetization near the interface. We 
calculate the change in magnetization from the pristine 
to the 0-vac and 0-ads systems and to the 20H and 
2H systems, AM, as ±1.7 /i^ and ±1.5 /iB in the su- 
percell, equivalent to ±0.42 e/S and ±0.37 e/S respec- 
tively (the remaining 0.1 electrons/holes stay in BTO, 
see Appendix). This extra electron density (which cor- 
responds to the electric displacement, because of the 
half-metallic nature of LSMO) resides in the interface 
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FIG. 4: Top: Schematic illustration of the change in band 
offset, ^vo, with polarization reversal. Bottom: Tunnel elec- 
troresistance (TER) vs BTO thickness. Experimental points 
taken from Garcia et al} (squares with dashed line fit) are 
compared with a theoretical expression^^ which uses the tun- 
neling barrier height expected from the BTO bulk polariza- 
tion, Po (solid line). Inset: Calculated band offset against 
electric displacement field for the three BTO states. The 
straight line fit is used to determine the band offset (and hence 
barrier height) at =bPo for the TER plot. 



region, decaying into the electrode with an associated 
Thomas-Fermi screening length (see Fig. 3). This situa- 
tion is similar to the carrier-mediated magnetoelectricity 
already predicted at SrTiOs/SrRuOa interfaceJ^ and in 
LSMO/BTO superlattices"^. In agreement with Ref.f27la 
competing interface antiferromagnetic type phase (called 
Al in Table 1 of Ref. |^ was found for the outwards 
BTO polarization. A similar magnetoelectric effect has 
recently been experimentally realizecP^^. 



E. Tunneling Electroresistance 

We now discuss how the electrochemical switching 
process may lead to the giant TER observed in the 
LSMO/BTO system^. In the simplest semiclassical ap- 
proximation, TER has an exponential dependence on 
the tunneling barrier shape^. The interface dipole, 
and hence band offset (^vo = ^vbm — ^f), at a 
metal/ferroelectric interface depends linearly on the elec- 
tric displacement field, 1), in a way that can be expressed 
with an effective screening length^^ Ag//. For LSMO- 
BTO we calculate Ag// = 0.11 A. Using the calculated 
values of the band offset (Fig. 4 inset) and the experi- 
mental band gap of BTO, we obtain the change in barrier 



height upon complete polarization reversal (for D = ±Po 
the potential in BTO is flat, i.e. the tunneling barrier 
shape is rectangular), Acp^ and the average barrier height, 
^ = (^out + ^in)/2. These values then yield an estimate 
of the TER using the exponential dependence^ on the 
barrier thickness, for large TER, 



TER; 



; exp 



(1) 



Fig. 4 compares this estimate with the experimental 
datsP showing that this simple model captures remark- 
ably well the essential physics of TER in this system. We 
note a recent study reported comparable shifts in Eyo 
(measured using photoelectron spectroscopy) on a similar 
ferroelectric/LSMO system upon polarization reversaP^. 
The origin of electroresistance effects in oxide nanotubes 
has also recently been suggested as redox reactionJ^. 
However the redox arguments there are fundamentally 
different - it is proposed that the electrons yielded by 
oxygen vacancies are directly available for conduction. 



IV. CONCLUSIONS 

In conclusion we have studied an electrochemical mech- 
anism for ferroelectric switching in thin films and pro- 
posed it as the origin of switchable ferroelectricity, TER 
and magnetoelectricity in a prototypical system. This 
work opens several avenues for future research. From 
the experimental point of view, it would be interesting 
to investigate the composition of a ferroelectric surface 
before and after switching (e.g. via the AFM tip), to 
verify whether reduced or oxidized gas-phase species are 
present (as suggested by our results). Also, this point 
could be indirectly checked by performing the AFM- 
mediated switching experiments in a controlled atmo- 
sphere, in analogy to the experiments of Bi et al. [19] 
on LAO/STO. From the theoretical point of view, a 
natural next step would be to perform a more detailed 
thermodynamic analysis of the stability of a ferroelectric 
surface (either pristine or decorated with adsorbates). 
This would involve exploring different coverages, possi- 
ble inhomogeneous polarization states, and the effect of 
temperature and other external perturbations. We hope 
that our results will stimulate further investigations along 
these (and possibly other) directions. 
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FIG. 5: (Color online) Spin-resolved layer- by- layer density of states centered around Ef for the LSMO/BTO systems. Positive 
DOS represents majority spin, negative DOS minority spin. Only the BTO layers (and 1 LSMO layer) are shown for clarity. 
Panels (a)-(e) correspond to pristine, 0-ads, 0-vac, 20H and 2H systems respectively. Panel (f) shows a profile of pfree through 
the various LSMO/BTO systems. 



Appendix 

Here we provide details of the electronic structure of 
the various LSMO/BTO systems. Figure 5 (a)-(e) shows 
the spin-resolved layer- by- layer density of states for the 
pristine, O-ads, O-vac, 20H and 2H systems. As dis- 
cussed in the main text, the electric displacement and 
polarization within BTO, D and P, and hence the va- 
lence band offset, Eyo, depend only on Q/S, the sur- 
face defect charge density, and not the surface chemistry. 
However, as discussed in Ref. 35, this is not strictly the 
case once Eyo becomes negative or reaches the band gap 
of BTO. At this point electrons or holes "spill out" in to 
the BTO layer. This "charge spill out" regime is favored 
by DFT, which often underestimates the band gap, and 



therefore can be an artifact of the calculation. 

In our case, due to the presence of a free surface, there 
is a further issue that was not explicitly considered in 
Ref. 35, i.e. the effect of surface states. In many cases, 
these fall within the bulk band gap of the ferroelectric 
film, and might cross the Fermi level of the metal, thus 
causing a significant spill-out of charge even when the 
bulk electronic bands are not directly affected. Note that 
surface states in ferroelectrics typically have a marked 
localized orbital character (either the atomic orbital of 
an adsorbate, or the 3d orbit als of the transition metal 
cation). Therefore, it is reasonable to suspect that DFT 
might introduce systematic errors in their ionization en- 
ergies (similarly to the energy location of the bulk band 
edges discussed in Ref. 35), and the metallization of a 



FIG. 6: Relaxed surface LSMO/BTO structure. Top: view along the [100] direction. Bottom: birds-eye view along the [001] 
direction. Only the top BTO layer is shown for clarity. Sr (blue), Ti (cyan), O (red), 0-ads (orange). Panels correspond to 
0-ads, 0-vac, 20H and 2H systems respectively (left to right). 



surface state should be regarded with analogous caution 
(for a detailed discussion of charge transfers at surfaces 
see, e.g., Ref. l36]) . 

Refs. I35l and [36l prescribe an analysis of the hole- and 
electron-like charge spill out. Following this prescrip- 
tion, we determine the free electron density profile, pfreei 
within BTO in Fig. 5(f) (using Eq. 25 and 26 of Ref.l35| 
and Eq. 19 of Ref. 36 for hole spill out). Out of all 
the five systems, the pristine one is unaffected, the neg- 
atively polarized (0-vac and 2H) systems are affected by 
electron spill out into the conduction band, and the two 
positively polarized (0-ads and 2 OH) ones are affected 
by hole spill out into surface states (see Fig. 5). In both 
0-vac and 2H pfree amounts to approximately 0.03 elec- 
trons per unit cell of BTO , which is a fairly mild effect 
(compare with approximately 0.15 electrons per unit cell 
in the KNO/SRO system of Ref. 35). In the case of O- 
ads and 2 OH, the surface 0{2p) states accommodate a 
total of approximately 0.1 holes per surface perovskite 
unit (Fig. 5(f)). Of course, estimating to what degree 
this charge spill is problematic, depends not only on the 
magnitude but also the purpose of the calculation. The 
charge spill out induces an error in two quantities that 
are discussed in this work: the total injected charge into 
the LSMO electrode, and the band alignment. Consider- 



ing the total injected charge, the impact of this error is 
trivial to estimate. In fact, 0.03 electrons per cell times 
N ^ number of BTO cells, corresponds exactly to the dif- 
ference between the actual induced spin in LSMO and 
the "ideal" limit of 2 Bohr magnetons per cell. This ob- 
servation can be directly used to estimate the error in 
the calculated band offset. In fact, we can assume in a 
first approximation that the band offset is linear in the 
electric displacement of the BTO cell adjacent to the in- 
terface, D inter- Usiug the above numbers for the 0-vac 
system, 

Dinter = 2e / S - 0.036(^-/5) = M/S = 1.7e/5, (A.l) 

where M is the induced spin, and S is the supercell sur- 
face area (or reciprocal of the defect density). This pro- 
vides an accurate estimate of the actual electric displace- 
ment "felt" by LSMO. Using this information, therefore, 
we can make a very accurate estimation of the linear band 
offset dependence with D. We used this analysis to make 
the plot shown in Fig. 4 inset of the main text. There- 
fore, whilst we understand the limitations of DFT, in this 
case they do not affect significantly our conclusions. 

Figure 6 displays the relaxed atomic structures of the 
BTO surface with 0-ads, 0-vac, 20H and 2H. 
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